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We consider a closed Friedmann-Robertson-Walker Universe driven by the back reaction from 
a massless, non-conformally coupled quantum scalar field. We show that the back-reaction of the 
quantum field is able to drive the cosmological scale factor over the barrier of the classical potential 
so that if the universe starts near zero scale factor (initial singularity) it can make the transition to 
an exponentially expanding de Sitter phase, with a probability comparable to that from quantum 
tunneling processes. The emphasis throughout is on the stochastic nature of back reaction, which 
comes from the quantum fluctuations of the fundamental fields. 

I. INTRODUCTION 

In this talk, based on work done in collaboration with Enric Verdaguer Q, we shall discuss a model of the Early 
Universe where a spatially closed, Friedmann - Robertson - Walker Universe avoids recollapse and launchs into 
inflationary expansion due to the effects of the back reaction of quantized matter fields. Central to the argument is 
the fact that this back reaction is both memory dependent and to some extent random, due to quantum fluctuations 
of the fundamental fields. The stress in particle creation and noise is the main novelty of our approach with respect 
to the by now large corpus of semiclassical cosmology || . 

To the best of our knowledge, this work is the most ellaborate application to date of the basic idea that quantum 
fluctuations of fundamental fields act on the geometry of the Universe as a stochastic energy -momentum tensor, put 
forward by several researchers |||-||]. We shall pass briefly on the technical details, which are contained in the original 
work, to concentrate on the physical ideas. In next Section, we give a general discussion of why and how stochastic 
terms ought to be included in Einstein equations; the following Section describes the model and the development of 
its solution, and we conclude with some brief Final Remarks. 



II. THE SEMICLASSICAL APPROXIMATION: PARTICLE CREATION AND NOISE 



Ever since the development of quantum mechanics and relativity theory early in this century, their final unification 
has been one of the most sought after prizes of theoretical physics. Moreover, as Hubble's observations first, and then 
the discovery of the cosmic microwave radiation, has taught us that the structure of our Universe is determined by 
events at the very beginning of its evolution, this search for unification acquired more than academic interest, since 
these early stages presumably demand a quantum description. However, in spite of some progress, specially in the 
sixties, the goal of a quantum theory of gravity seems now as elusive as ever p0| . 

In the meantime, it has been realized that even without a full theory it was possible to find answers to most 
pressing questions, at least from the cosmological and astrophysical viewpoints. Basically the theoretical expectation 
is that quantum gravitational effects will become relevant at the Planck scale (of 10 19 GeV in natural units), while 
quantum effects associated to matter are conspicuous at the scales associated to the masses of those particles, orders 
of magnitude below that. Therefore, it makes sense to develop models where the geometry of the Universe (or, say, 
the space surrounding a compact object) is treated in the terms of Einstein's theory (which is to say, classically) while 
matter is quantized ||. As shown by Parker jllj], even these restricted models leave ample room for exotic behavior, 
probably its most spectacular manifestation being Hawking radiation [T^jr^ |. 

Now, there being no question that quantum fluctuations weigh, as shown by the Casimir effect, its presence should 
affect the evolution of the said geometry, on equal footing with other forms of matter, gravity included. So, although 
models of quantum test fields on given backgrounds have been immensely useful to clarify the formal aspects of the 
theory, the real goal is to develop self consistent models, where the feedback between the quantum matter and the 
classical geometry is fully accounted for. This raises the issue of how the quantum matter affects the geometry. 

The simplest and earliest answer to this question is that it is possible that while the field amplitudes associated to 
matter are to be described by operators in some functional space, other observables such as energy - momentum, which 
are composites of field operators, may develop a condensate or c-number part. In the simplest case, this condensate 
may be computed as the expectation value of the corresponding composite operator in an adequate quantum state, 
and this expectation value is to be included as the quantum matter contribution to Einstein's equations p3 . 
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By now, there is a substantial body of work showing that this approach can be made mathematically consistent 
. However, there are also serious doubts concerning whether it is physically satisfactory Jl^,|l6|. One can formulate 
these doubts with any degree of sofistication, but the basic argument goes as follows. 

In a typical self - consistent evolution problem, the identification of the energy - momentum tensor of the matter 
fields involves a normal ordering or subtraction procedure, to take care of the ever present divergences of quantum 
field theory. Suppose this normal ordering amounts to substract the energy - momentum defined with respect to some 
local vacuum state, the most sensible choice being a local adiabatic vacuum of some order. Suppose we let the model 
evolve from some time t, to some final time tf. In almost every interesting case, the evolving geometry will mix the 
positive and negative frequency components of the field operator, so that the creation and destruction operators of 
the adiabatic model at time tf will be related to that at time U by a Bogolubov transformation. This transformation 
is characterized by two complex parameters, a and /3, with \a\ 2 — |/3| 2 = 1. If at time t, a given one particle state was 
occupied by n adiabatic particles, then at time tf we find, in average, n' = |/3| + n \a\ particles there (we assume 
Bose statistics, for concreteness). This average is reflected in the evolution of the mean energy -momentum. But the 
state at times tf will not, in general, be a state with a well defined occupation number (certainly not, if the state at 
time ti was), and when we look at the dispersion in particle number, we see that unless |/3| 2 -C n, in which case we 
did not need to bother with self consistency, the average fluctuation in particle number is of the order of the mean 
occupation number itself. So the average energy momentum of the created particles is only vaguely related to what 
might have happened " on the spot" . 

The issue then arises of how to introduce these fluctuations in an actual model. Maybe the final answer is that 
nothing short of full quantum gravity is truly satisfactory, but if one wishes to retain the classical character of the 
metric, then it seems that the only possibility is to add, to the mean energy - momentum of the fields, a stochastic 
component, which would represent the leading effects of the quantum fluctuations. 

J. Halliwell and others fl7],|l8| have developed a picture which shows that this classical, stochastic energy momentum 
at least makes sense. In this approach, the classical geometry is seen as an apparatus which measures the energy - 
momentum of the quantum matter fields, and reacts to the measured value. It is then seen that the results are c - 
numbers, but they progression in time does not follow deterministic laws. It is natural to extract the deterministic, 
mean evolution, and call the remainder random noise. In the limit where the measurements are inaccurate enough, 
and repeated often enough, we may assume that the evolution of the fluctuations will not be greatly perturbed by 
the measurement process, and so they may be computed from ordinary quantum field theoretical rules. 



A. The closed time path effective action 

In this limit, therefore, a new paradigm appears, derived from the work of Feynman and Vernon E3] . We now 
regard the geometry as an open system evolving in the environment provided by the matter quantum fluctuations 
[j20|. Since the detailed evolution of the later is deemed irrelevant, our only concern is to estimate the influence 
action, namely, the modification to the gravitational action due to the influence of the environment. In the limit 
where the geometry is actually taken to be classical, the task becomes identical to that of computing the so-called 
Schwinger - Keldysh effective action |2lJ . We shall disregard the somewhat technical distinctions between these two 
objects, regarding them as the same; in a nutshell, the Schwinger - Keldish or closed time path (CTP) effective action 
(EA) is the influence action evaluated over an infinite time range, which actually takes care of some difficulties in the 
evaluation of the influence action over finite lapses. 

The CTPEA is a truly remarkable object, which achieves the miracle of providing a well defined variational method 
to derive causal but non local in time equations of motion. It is not hard to see where the difficulty lies. Suppose 
you have a system described by some variable (ft if) > and write for it an action functional S [(ft] whose variation yields 
the equations of motion SS/6<p (t) — 0. If the equations are causal, then S 2 S/5(j) (t) 5(ft (t 1 ) = whenever t' > t. Since 
second derivatives commute, the second derivative actually vanishes if only t' ^ t, and the action must be necessarily 
local in time. Since all dissipative effects are physically limiting cases of non local time interactions (when the response 
time of the bath is much shorter than the characteristic time of the system), it follows more generally that there is 
no variational principle for dissipative, causal evolutions. 

The CTPEA achieves the impossible by adding to every degree of freedom </> + (t) a mirror degree of freedom (ftr (t) , 
so that the CTPEA T = T[<ft + ,(ft~] and the equations of motion are ST/S(ft + = 0. The right count of degrees of 
freedom is restored by imposing, after the variation has been taken, the constraint <fi + = (ftr = (ft, the physical degree 
of freedom. Causality only demands 

when t < t' (1) 



5(ft+ (t) 5(ft+ (t>) 5(ft+ (t) deft- (f) 
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For example, suppose that the solution to the equations of motion is just <fi = 0, and we seek the dynamics of small 
fluctuations. The only quadratic action compatible with the causality constraint has the form 

T = \J dtdf {M(^(M')W(i') + M(t)^(M')M(i')} (2) 

where [(f)] = cj) + — <fi~; {4>} = cf) + + <j>~ , and T> (t, t') = if t < t' , The equations of motion are 

d£ V (t, t') (V) = (3) 

and we see that there is no obstacle to causality, with no further restrictions on locality. 

To actually compute the CTPEA, Schwinger observed that the mean value of the field could be obtained from a 
generating functional 

Z[J+,J-} =j w l J+ > J -]=(f(e- i f J ~*)T(e i f J+ *)) (4) 

where $ is the Heisenberg field operator, the expectation value is taken with respect to the corresponding quantum 
state, and T, T stand for temporal and antitemporal ordering, respectively. Indeed, if we define 

Then in the limit J ± — » 0, we obtain (f> + = <f)~ = (<&) . This suggests defining F as the Legendre transform of W 

-]=W[J+,J-]- J (J+0+-J-0-) (6) 



So in general we obtain 



and for the physical mean field 



sr 



— =TJ± 



= (7) 



as required. Observe that on top of Eq. (]]]), the quantum CTPEA obeys T [4> + , <f)~] = — T [<f)~ , (j> + ]* , so that the 
kernel T> must be real (and so are the equations of motion) and Af is pure imaginary (J\f — iN) . 

The CTPEA admits a functional representation derived from the usual background field methods p3] 



(8) 



The variables of integration must coincide at some very large time in the future, and we have included the 
information on the quantum state in the integration measure. In our case, we have both system (ip) and environment 
(ip) degrees of freedom, but we associate a mean field only to the former. The classical action S will be the sum of 
the system action S s , the environment action S e , and the interaction term Si. In the semiclassical limit, we neglect 
the deviation from the mean of the system variables, and the expression for the CTPEA simplifies to 

e v[* + ,*-] =e i[s.(*+)-s.(*-)] J D ^ Dr expz[5 e (^)+5 4 (0+>+)-5 e (^-)-5. i (r,^)] (9) 

where, again, the integration is nontrivial thanks to the integration measure and the future boundary conditions. 

Eq. provides a well defined recipe for the CTPEA, and from thence finding the effective equations of motion 
is only a matter of computing power. The question arises of when the solutions to the semiclassical equations are 
relevant to the description of our experience. Mostly what we wish to do is to compute expectation values of system 
observables. In so far as these observables do not involve environmental variables, their expectation values will admit 
representations such as 
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A = J D<j)+D<i)- A [<P+, <j)-] e lI > <* \ (10) 

which under the saddle point approximation reduces to 

A = A[4>,<f>] (11) 

4> being the solution to the mean field equations. In other words, the semiclassical equations will be useful if our 
only interest is to compute expectation values for system observables, such as the integral Eq. ([ll]) may be evaluated 
by saddle point methods. 

The CTPEA is a powerful method to derive equations of motion for open (or effectively open) systems, which are 
guaranteed to be both real and causal [^3| , ^4| . It is possibly the most powerful method at hand to study nonequilibrium 
evolution of quantum fields, specially if combined with more sophisticated resummation methods, which allow us to 



keep track of higher Schwinger functions along with the mean field 1 25 26|] . But two things ought to send an alarm 
signal. First, if we only care about the equations of motion, we are only using a small part of the information encoded 
in the CTPEA; for example, the kernel M is irrelevant to the linearized equations (||). Second, there is no noise in 
the mean field equations. We must stress that this is not only unseeming, but, insofar as the equations of motion are 
dissipative, it is actually wrong, as it violates the necessary balance between fluctuations and dissipation fl27j. 



B. Where is the noise? 



To the best of our knowledge, it was Feynman who pointed out that the untapped terms in the CTPEA contained 
the information about noise jl9|]. For simplicity, assume the CTPEA has the quadratic form Eq.(Q). Then we have 
the identity 



ir [<^1 = J Djp\j]e X ph J dtdt' [4>] (t)v(t,t/){<f>}(t') + i Jdtji 



(i) (t) (12) 



with 

(j (t) j (t')) =jDjp [j] j (t) j (*') = N (t, f ) (13) 
Now we can write the average Eq. ( ^fj| ) 

A= J Dj p [j] J D4> + D4>- A [<f>+, 4T] exp || j dtdt 1 [</>] (t) V (t, t') {4>} (*') + J dt j (t) [cp] (t) j (14) 
And use saddle point evaluation in the inner integral 

A = J Dj p\j]A[^^ ] (15) 

where 4>j is the solution to 

' dt'V(t,t')<p(t') = -j(t) (16) 



In this way, we have transformed the average of the observable A over the quantum fluctuations of the environment 
into the ensemble average over the realizations of the "noise" j(t), at the same time upgrading the semiclassical 
equations to the Langevin equations Eq. (|l^). As expected, the relevant information on the noise (its correlation 
function), is given by the "useless" part of the CTPEA, namely the kernel N = — iV. 

Of course, eq. ( ]l2] ) is not the only way to decompose the CTPEA into partial integrations. The point is that 
this particular decomposition makes physical sense. To see this, assume the interaction term in the action takes the 
particular form 



dt S <p (17) 
with (S) = when = but otherwise arbitrary. Then |0] 
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N(t,l/) = ~({E(t),Z(t')}) (18) 



the expectation value being computed at 4> = 0. Indeed, the Heisenberg equation for this model is 

6S» 



(19) 



We assume that the Heisenberg operator for the system variable is close to a c-number. Also, in the presence of a 
non zero background </>, the operator S will generally develop a nonzero expectation value (S)^. Subtracting this, we 
get 



(S-(H)J (20) 



The CTPEA, if we forget about noise, leads to the same equation with no right hand side. As discussed above, this 
is not acceptable. So the question is, what is the sensible way of replacing the q- number operator in the right hand 



side of Eq. (20) by a c-number stochastic source. Of course, some loss of information (specially concerning quantum 
coherence) is unavoidable, and in particular cases this may invalidate the whole procedure. But whenever quantum 
coherence is not the main concern, a Gaussian source with self correlation as in Eq. ( |l8| ) (we again neglect (S)^ , 
which vanishes at (f> = 0, since we assume small fluctuations, but this is not essential) is the time honored answer, 
and indeed the only answer compatible with the fluctuation dissipation theorem p8| . 



III. NOISE INDUCED INFLATION: WHEN NOISE MATTERS 



As we have seen in the previous Section, the CTPEA provides a systematic framework in which to study semiclassical 
evolution, taking into account at least the leading effects due to quantum fluctuations of matter fields. However, 
putting this framework into work is by no means a simple task. Indeed, after the observation that the CTPEA 
provides a simple way of deriving Einstein - Langevin equations ||, it was realized that actually carrying out the 
derivation is a research project in itself |§-f§, to say nothing of solving those equations once derived ||. So it is 
natural to wonder if noise makes such a difference as to justify this trouble. 

The basic problem is that, while it would be easy to find problems where the noise level is huge, this same noisyness 
would lead to the suspicion that the whole semiclassical approach is breaking down. The real challenge is to find a 
problem where the semiclassical approximation is reliable, and still noise makes a difference. Indeed, in the case of 
conventional, noiseless semiclassical theory, such a problem is Hawking evaporation of large black holes: a weak effect, 
which puts in no jeopardy the validity of the semiclassical approximation, but whose result is utterly impossible in 
terms of the classical theory alone. 

In this talk, I will report on one such problem, a cyclic Universe, provided with a cosmological constant but prevented 
from inflating by the potential barrier from its own spatial curvature. In each cycle, the semiclassical effects induce a 
transition form one classical orbit to another; the change is small for each cycle, but overall it offers an escape route 
with no classical analog. Our goal is to compute the average escape probability due to semiclassical effects. 

This problem has an important precedent, the calculation of the tunneling amplitude due to Vilenkin ps|| . It is 
important to realize the similarities and differences between these two approaches. Vilenkin's calculation was fully 
quantum gravitational, but it contemplated only the effects of the gravitational field. It was tacitly assumed that, 
if any matter fields were present, they would at most affect the prefactor of the exponentially suppressed tunneling 
probability |50] ^||. Our calculation is only semiclassical, but we put the stress precisely on the effects of the matter 
fields. From the point of view of the usual instanton approach, we could say ours is a highly nonperturbative 
evaluation of the prefactor, since we go well beyond the test field - one loop approximation. The result is that the 
escape probability due to the fluctuations in matter fields is at least as large as the tunneling one, suggesting that 
in Nature both must be taken into account. We shall not discuss subsequent developments related to Vilenkin 's 
proposal pq . 

Since the calculation of the escape probability due to semiclassical effects is discussed in some detail elsewhere [|l] , 
here I shall only give a general discussion of the several steps involved, the peculiar difficulties of each one, and how 
they could be overcome. 
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A. The model 



Our model is based on a spatially closed, homogeneous Friedmann - Robertson - Walker (FRW) model, with a 
metric 



ds 2 = a 2 (t) [-dt 2 + g ij {x k )dx i dx i ) , k = 1, . 



1- 



(21) 



where a(t) is the cosmological scale factor, t is the conformal time, and gu{x r) is the metric of an (n — l)-sphere of 
unit radius. Since we will use dimensional regularization we work, for the time being, in n-dimensions. Matter is 
described by a quantum scalar field $(x^), where the Greek indices run from to n — 1. The classical action for this 
scalar field in the spacetime background described by the above metric is 

n- 2 



S m = - / dx n ^g 



4(n- 1) 



(22) 



where <?oo — fl2 j 9oi — 0; 3ij — a2 9ij> 9 is the metric determinant, v is a dimensionless parameter coupling the field to 
the spacetime curvature (y — corresponds to conformal coupling), R is the curvature scalar which is given by 

R = 2(n- 1)4 + (n - - 4) 4 + (« " " 2 )4> ( 23 ) 
a* a* a z 

where an over dot means derivative with respect to conformal time t. Let us now introduce a conformally related 
field * 



the time dependent function U (t) 



U(t) = 



2 (t)R(t), 



and the d'Alambcrtian □ = —d 2 + A'" x ^ of the static metric ds 2 = a 2 ds 2 . The action may be written as, 



5/1 



dtdx 1 . . . dx r ' 



(n-2) 2 



(24) 
(25) 

(26) 



When v = this is the action of a scalar field f in a background of constant curvature. The quantization of this 
field in that background is trivial in the sense that a unique natural vacuum may be introduced, the "in" and "out" 
vacuum coincide and there is no particle creation ^|. This vacuum is, of course, conformally related to the physical 
vacuum, see (^4|). The time dependent function U(t) will be considered as an interaction term and will be treated 
perturbativelly. Thus we will make perturbation theory with the parameter v which we will assume small. 

To carry on the quantization we will proceed by mode separation expanding ^(x 11 ) in terms of the (n— l)-dimensional 
spherical harmonics Yl(x l ). The coefficients ^c(i) are just functions of t (1-dimensional fields), and for each set (I, k) 
we may introduce two real functions ^g(i) and <pj-{t) defined by 



(27) 



$L = — ( 6l 

k ~ V2 V k 

The action becomes the sum of the actions of two independent sets formed by an infinite collection of decoupled 



time dependent harmonic oscillators 

1 

2 



s - = \ /*EE f(4) 2 - M i (4) 2 + m 
k=i t l 



where the dots stand for an identical action for the real 1-dimensional fields <?v(i)- 

We will consider, from now on, the action for the 1-dimensional fields <fig only, 
dimensional fields yAt) are, from 



(28) 



The field equation for the 1- 



(29) 



which in accordance with our previous remarks will be solved perturbatively on U it) . The solutions of the unperturbed 
equation can be written as linear combinations of the normalized positive and negative frequency modes, fk and f£ 
respectively, where 

f k (t) = ^=exp(-iM k t). (30) 
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B. Closed time path effective action 



We are now in the position to compute the regularized semiclassical CTP effective action. This involves a careful 
consideration of the infinities arising in perturbation theory, but after the dust settles, the result is 

TcTpla*] = S« m [a+] - S^ m [a~] + S? F [ a ±], (31) 

where the regularized gravitational and classical matter actions are, 

Sfjfl] = J dtQa 2 + lj - 2tt 2 J dt^A* + dtUi(t)ln(a f i c ). (32) 

and the influence action 

5f i ,[a ± ] = i J dtdt' AU(t)H(t - t'){U{t')} + % -J dtdt 1 AU(t)N(t - t')AU(t'), (33) 
where we have defined 

AU = U + -U~, {U} = U + + lf-. (34) 

Computing the kernels H and N involves the consideration of Feynman graphs, where the internal legs represent 
propagators for a particle in a closed space. This calculation may be carried out exactly, but the result is that, 
unless for orbits with very small amplitude, the effect of spatial curvature is not really important. It is convenient 
to compute these kernels as in a spatially flat FRW Universe with the same radius, which amounts to consider a 
continuous, rather than a discrete, spectrum of modes. The result is 

f°° 7T 

N(u)= dkcos2ku= — S(u). (35) 
Jo 16 



H(u) = ~Pf 



9(u) 



H o(u). 



(36) 



The distribution Pi(0(u)/u) should be understood as follows. Let f(u) be an arbitrary tempered function, then 

-0(u) 



duPf 



/(«) 



lim 



du^ + /(Q)lne 



(37) 



The approximation of substituting the exact kernels by their flat space counterparts is clearly justified when the 
radius of the universe is large, which is when the semiclassical approximation works best. 

. The imaginary part of the influence action is known |^[^]3"7j-{4C| to give the effect of a stochastic force on the 
system, and we can introduce an improved semiclassical effective action, 



1 



S eff [a ± ;Z}=S* [a + }-S* m [a-] + - / dtdt' AU(t)H(t - t'){ U(t')} + / A 



where is a Gaussian stochastic field defined by the following statistical averages 

«(*))= o, mm) = N(t-t>). 



(38) 



(39) 



The kernel H in the effective action gives a non local effect (due to particle creation), whereas the source £ gives 
the reaction of the environment into the system in terms of a stochastic force. 
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C. The Einstein - Langevin equation 



The dynamical equation for the scale factor a(t) can now be found from the effective action (|3g) in the usual way, 
that is by functional derivation with respect to a + (t) and then equating a + = a~ = a. These equations include the 
back-reaction of the quantum field on the scale factor; they improve the semiclassical equation by taking into account 
the fluctuations of the stress-energy tensor of the quantum field [ p^|Jl^ ] . However, they also lead to the typical non 
physical runaway solutions due to the higher order time derivatives involved in the quantum correction terms. 

To avoid such spurious solutions we use the method of order reduction plj j. In this method one assumes that the 
equations obtained from the CTPEA are perturbative, the perturbations being the quantum corrections. To leading 
order the equation reduces to the classical equation, which, in terms of scaled variables 



m^a(t), A.^A*. (40) 

reads 

b + b (l - ^A6 2 j = 0(v). (41) 

The terms with b or with higher time derivatives in the quantum corrections are then substituted using recurrently 
the classical equation (Ell). In this form the solutions to the semiclassical equations are also perturbations of the 



classical solutions. Thus, by functional derivation of (38), we can write the stochastic semiclassical back-reaction 
equation as 

p = -V (b) - 5V'{b) + F(b,p,t) + J(S,b,p), (42) 
where a prime means a derivative with respect to b, and we have introduced p = b. The classical potential V(b) is 

V{b) = \b 2 -^b\ (43) 
An schematic plot of this potential is given in Fig. 1. 

The remaining terms in eq. (E3) represent the quantum corrections. The first one is purely local 



-b 2 - —b 4 - p 2 ]n(bn) 
2 48 y K w 



(44) 



where we have already implemented order reduction. 

From this point on, we shall disregard the local quantum correction to the potential, SV(b). In the region where 
semiclassical theory is reliable, this is only a very small correction to the classical potential; moreover, we are concerned 
with such phenomena where the semiclassical behavior is qualitatively different from the classical one, which is not 
the case for these corrections. 

The term F(b,p, t) involves nonlocal contributions and may be written as, 

dU d 2 fdU \ f d 2 1 b ) 

F(b lP ,t) = -^1-0 = *»{ W - b - V 2 ) 7 < ^ 

where I(b,p,t) is defined by 

/oo 
dt'H(t - t')U{t'). (46) 
-oo 

and 

U{t) = -&u\"- + l\. (47! 




After order reduction, U (f) must be evaluated on the classical orbit with Cauchy data b (t) — b, p (t) — p, whereby 
it reduces to U = — Kvb 2 . Observe that, in fact, this approximation makes the equation of motion local in time, 
though non longer Hamiltonian. Finally, the function J is the noise given by 
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and, after order reduction, by 



| 2jp 2ZV'(b) 2t;p 2 
b b 2 b 2 b 3 



(48) 



with £ (t) denned in (p9|) in terms of the noise kernel. 



D. The classical orbits 

Before continuing, it is convenient to pause and consider the classical orbits, as described by Eq. ([id]). They 
represent a particle moving in the one dimensional potential well Eq. (p3[), plotted schematically in Fig. 1. This 
evolution preserves the Wheeler - DeWitt operator (energy, for short) 

H(b,p) = ±p 2 + V(b) (49) 

Physically, the value of H on a given orbit is the energy density of radiation present besides the cosmological 
constant. 

Fig. 2 is a schematic representation of the classical phase space. There is a stable fixed point (nothing) at b = p = 0. 
This is the starting point of Vilenkin's calculation. There is an unstable fixed point corresponding to 

H = E s = ^ R -b = 2^W s -p = Q (50) 

This is an Einstein type static Universe. When H > E s , orbits are free to expand forever. 

For H < E s , we have two types of orbits. Those outside the well contract at first, until they reach the classical 
turning point b+ and bounce off. The Dc Sitter Universe, with H — 0, belongs to this family. In their final stages, 
these orbits are essentiallly identical to the ever expanding ones (cosmic no hair theorem). Orbits inside the well with 
< H < E s bounce ethernally between the turning points ±6_ (there is no problem with a negative radius of the 
Universe, since only b 2 has a physical meaning; we may also think of b = as a perfectly reflecting boundary). The 
actual location of the turning points is 



6± = ±E S 



(51) 



The frequency of oscillation is 1 for H <C E s , and vanishes as H — > E s . This limiting value corresponds to two 
orbits, exponentially departing from and approaching to the unstable fixed points, the so-called separatrices. 

To fix ideas, let us adopt for A a value consistent with Grand Unified scale inflation, which in natural units means 
A ~ 10 -12 . Then the value of the separatrix energy is very high, E s ~ 10 12 . Our problem is to find a way for a typical 
Universe (H ~ 1), trapped within the well, to climb out of it and inflate. As we shall see, this is possible thanks to the 
combination of the diffussive effect of quantum fluctuations and a runaway instability associated to particle creation. 

It ought to be clear that the final state of this evolution will be very different than in the instanton approach. In 
the quantum claculation, the Universe emerges from under the barrier as an empty, H = 0, De Sitter Universe. In our 
calculation, the Universe goes above the barrier, and emerges with a large amount of radiation H — E s , corresponding 
to particles created while inside the well. Physically however the difference is minor, as this energy gets diluted in a 
few e-foldings by the inflationary expansion. 



E. From Langevin to Kramers 



Now we want to determine the probability that a universe starting at the potential well goes over the potential 
barrier into the inflationary stage. The magic of the CTPEA has turned an originally quantum problem into a 
statistical mechanical one, indeed a classic problem associated to the name of Kramers |42]. Observe that we are not 
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interested in the features of solutions associated to peculiar realizations of the noise, but rather on a noise averaged 
observable. Therefore, it is convenient to perform the noise average at the outset, introducing the distribution function 



f(b,p,t) = (6(b(t)-b)S(p(t)-p)), (52) 

where b(t) and p(t) are solutions of equation (|42j) for a given realization of £(t), b and p are points in the phase 
space, and the average is taken both with respect to the initial conditions and to the history of the noise. After some 
standard manipulations we arrive at the so-called Kramers' equation |Jf3f 

where the curly brackets are Poisson brackets, i.e. 

{HJ} = -p(df/db) + V'(b)(df/dp), 

and 

This term will be called the diffusion term since it depends on the stochastic field 

We notice that in the absense of a cosmological constant, we get no diffusion. This makes sense, because in that 
case the classical trajectories describe a radiation filled universe. Such universe would have no scalar curvature, and 
so it should be insensitive to the value of v as well. 



F. From Kramers to Fokker-Planck 

In the usual statement of Kramers' problem, the system is described by a single variable x and obeys a Fokker-Planck 



equation 44 



df = d_ 

dt dx 



dx dx 



(55) 



where T is the temperature (which fixes the sign of the diffussion term) and F is the free energy (rather than a 
potential). Activation is studied from the properties of the steady solutions of this equation, and the answer is the 
so-called Arrhenius formula 

p „ £ -F^/T (56) 

where F max is the value at the peak of the free energy barrier. 

Our Kramers equation is certainly more involved than Eq. (^5|), because it describes other phenomena besides 
tunneling. Basically, there are three things going on. Given a generic distribution function /, its dynamics consists 
mostly on the representative phase space points being dragged along the classical orbits, with a time scale of the order 
of a typical period. On a larger time scale, we have the diffussion process, which makes / evolve towards a quasi 
equilibrium, steady solution. Finally, there is activation, on an even larger time scale. 

Since our concern is this third process only, it is convenient to get rid of the two faster ones. We get rid of classical 
transport by defining a new distribution function which counts the number of Universes on a given classical orbit, 
rather than on a phase space cell. This new distribution function does not tell us where in the orbit we are, but we 
do not need that to study activation. We achieve this by transforming the problem to action - angle variables, and 
averaging over the angles |^5| . Finally, we get rid of the approach to quasi equilibrium by assuming a steady solution 
from the beginning. 

The averaged Kramers equation becomes, 

o f _^a fu(j )d j s/ j 



dt 4 dJ [ fl dJ 

where 
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DU) = ^ / n d6b 2 p 2 , ( r,si 

s < J » = i?l Hs^'M (59) 

This equation may be written as a continuity equation dtf + djK = 0, where the probability flux K may be identified 
directly from (p7|). We see that, as in Kramers' problem, stationary solutions with positive flux Kq should satisfy 

— dJ- S{J)f = -^Al Ko - (60) 

From now on it is more convenient to use the energy £ as a variable instead of J, where E = H(J). D and S 
individually behave as E 2 times a smooth function of E/E s , and their ratio is relatively slowly varying. At low energy, 
we find D ^E 2 /2 and S ^E 2 /4. As we approach the separatrix, D — >0.96E 2 and S -+1.18E 2 . Meanwhile, the ratio 
of the two goes from 0.5 to 1.23. This means that we can write the equation for stationary distributions as 

where /3 and g are smooth order one functions. There is a fundamental difference with respect to Kramers' problem, 
namely the sign of the second term in the left hand side. In the cosmological problem, the effect of nonlocality is to 
favour diffussion rather than hindering it. We may understand this as arising from a feedback effect associated with 
particle creation (see pfl). 



G. The activation amplitude 

Fig. (3) is a out of scale, schematic plot of the solution of Eq. (|6l|). For E <C 1, the solution diverges as 1/E; for 
E > 1, and for twelve decades thereafter, it grows exponentially. Of course, our analysis doeas not hold beyond the 
separatrix, but it can be shown that / turns around there, decaying as a power of E as E — > oo . 

This behavior cannot be extrapolated all the way to zero as it would make / non integrable. However we must notice 
that neither our treatment (i.e., the neglect of logarithmic potential corrections) nor semiclassical theory generally is 
supposed to be valid arbitrarily close to the singularity. Thus we shall assume that the pathological behavior of Eq. 
( |6ij ) near the origin will be absent in a more complete theory, and apply it only from some lowest energy Eg ~ 1 on. 
There are still 12 orders of magnitude between Eg and E s . 

We may now estimate the flux by requesting that the total area below the distribution function should not exceed 
unity. Unless the lower cutoff Eg is very small (it ought to be exponentially small on E s to invalidate our argument) 
the integral is dominated by the peak around E s , and we obtain 

K < (prefactor) exp [-0(E a )E a ] . (62) 

The prefactor depends on A, v, g (1), j3 (1), a and the details of the peak shape, but only logarithmically on the cut 
-off. So we can take this as a bona fide prediction of noisy semiclassical theory. Using E s — 3/ (2A), (3 (E s ) = 1.23, we 
get 

Ko < (prefactor) exp (— — r— ) • (63) 



A 

This is the semiclassical result to be compared against the instanton calculations , which yield 

P~exp(-~). (64) 

We see that the probability of noise activation is indeed larger than that of quantum tunneling. 

A final pertinent question is, was noise truly necessary? After all, one could imagine there would be particle creation 
in each cycle, and the accumulation of these particles alone would make inflation progressively easier to achieve. Of 
course, in the absense of diffussion our argument would need to be entirely redrawn; for example, it is not longer 
clear than an arbitrary / would tend towards an steady solution, or on which timescales, nor that the steady solution 
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would have an acceptable behavior beyond the separatrix. However, one could try to see what happens if one simply 
kills the diffussion term in Eq. ( |6l| ) . The equation still admits a solution, but the divergence at the origin gets worse 
(1/E 2 ), and indeed the flux seems to depend solely on the value of the infrared cut - off, Kq ~ Eg, which cannot be 
predicted within the semiclassical theory. It is only the combination of particle creation and diffussion which sets up 
the mechanism by which we get a definite result. 

On the other hand, an infinity of cycles is not truly necessary. There is a similar result in models where the Universe 
is restricted to a single cycle; the activation probability is somewhat lower, but still higher than the tunneling one. 
We can understand this by analogy to the problem of black hole formation in a box, where a big hole can form by 
slow accretion of smaller holes, or by a sudden, single large fluctuation. The difference is, of course, that in the black 
hole case the latter is more likely than the former [j46| . 

IV. FINAL REMARKS 

We have reported on a cosmological process where quantum induced noise and particle creation combine to yield 
a behavior notoriously different than expected from classical theory, or even conventional, deterministic semiclassical 
gravity. The strenght of the effect is indeed comparable to a purely quantum calculation, which shows by the side 
that treating matter as test fields in quantum gravity may not be justified. We believe this work is meaningful in at 
least three different levels: 

a) of course, our results are most important as an step forward in the development of stochastic, semiclassical 
cosmological models. By now, the mathematical and physical basis of such models is rather well understood, but the 
development of actual models and the gathering of hard predictions is lagging behind. Our calculation has demanded 
the application of a number of techniques which are not common tools of the trade in cosmology, and could serve as 
an example for future projects. 

b) the relevance to cosmology may seem minor, since it does not seem likely that our Universe be spatially closed. 
However, the situation changes if the original question is rephrased as: is it possible that a horizon size, overdense 
region in the early Universe, with a homogeneous but subplanckian value of the inflaton field, may avoid collapse and 
inflate? Clasically, the answer is no, and this negative result may well be the bane of inflationary models |47|]. Our 
results suggest that semiclassically things may turn around. 

c) finally, it has been observed that noise and dissipation are generic to all effective theories (io) . So we must expect 
that similar results will be found in the analysis of nuclcation phenomena in other effective theories as well, specially 
in quantum field theories ]4S|| , and in out of equilibrium situations. Indeed, an approach such as ours seems to be the 
only way of analyzing tunneling in situations where the environment changes on timescales comparable to the time it 
takes to nucleate a bubble, an essentially virgin field right now. 

We continue our research on all these levels, and hope to report soon on new results. 
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V. FIGURE CAPTIONS: 



A. Fig 1. 

A schematic plot of the classical potential; it vanishes as b 2 when b — ► 0, it has a maximum, and decreases without 
bound for larger Universes. Classical evolution preserves the Wheeler - DeWitt operator H = p 2 /2 + V(b). If H 
exceeds the maximum of the potential, the corresponding orbit either expands forever or collapses to the singularity. 
For lower H, the classical orbit bounces off the outer classical turning point. For positive H below the maximum of 
the potential, we have periodic bounded orbits representing an eternal cyclical Universe. 

B. Fig. 2. 

A sketch of the classical phase space. Only half is shown, the other half being the mirror image. We can see the 
stable (elliptic) fixed point at the origin, and one of the unstable (hyperbolic) fixed points. The separatrices connect 
the unstable points to each other, and divide the region of periodical motion (within) from the region of unbound 
motion (outside). The normalization is b' = b/2\/E~ s , p' = p/\/2E s . 

C. Fig 3. 

A qualitative plot of the equilibrium distribution function, as a function of E/E s . We can appreciate the divergence 
towards the origin, the exponential rise towards the separatrix, and the falling off in the inflationary region. For generic 
values of the cutoff, the area under the curve is dominated by the peak at the separatrix. 
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